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ABSTRACT 

Context. A significant number of cataclysmic variables were detected as hard X-ray sources in the INTEGRAL survey, most of them 
of the magnetic intermediate polar type. 

Aims. We present a detailed X-ray broad-band study of two new sources, IGR J00234-I-6141 and IRXS J213344.1+510725, that allow 
us to classify them as secure members of the intermediate polar class. 

Methods. Timing and spectral analysis of IGR J00234+6 14 1 are based on a XMM-Newton observation and INTEGRAL publicly avail- 
able data. For IRXS J213344.1+510725 we use XMM-Newton and Suzaku observations at different epochs, as well as INTEGRAL 
publicly available data. 

Results. We determine a spin period of 561.64 ± 0.56 s for the white dwarf in IGR J00234+6141. The X-ray pulses are observed up 
to ~ 2 keV. From XMM-Newton and Suzaku observations of IRXS J213344.1-I-5 10725, we find a rotational period of 570.862 ± 0.034 
s. The observations span three epochs where the pulsation is observed to change at different energies both in amplitude and shape. In 
both objects, the spectral analysis spanned over a wide energy range, from 0.3 to 100 keV, shows the presence of multiple emission 
components absorbed by dense material. The X-ray spectrum of IGR J00234+6141 is consistent with a multi-temperature plasma 
with a maximum temperature of ~ 50 keV. In IRXS J213344.1+510725, multiple optically thin components are inferred, as well as 
an optically thick (blackbody) soft X-ray emission with a temperature of ~ 100 eV. This latter adds IRXS J213344. 1+510725 to the 
growing group of soft X-ray intermediate polars. 

Conclusions. These two intermediate polars, though showing similar rotational periods and being hard X-ray sources, appear differ- 
ent in several respects, among which the presence of a soft X-ray component in IRXS J213344.1+510725. This source also emits 
circularly polarized light in the optical band, thus joining as a fifth member the group of soft and polarized intermediate polars. How 
the hard X-ray and polarized emission are thermalized in these systems is a timely question. 

Key words. Stai-s: binaries: close - Stars: individual: IGR J00234+6141, IRXS J213344.1-I-510725 - Stai's: novae, cataclysmic vari- 
ables - X-rays: stars - Accretion, accretion disks 



1. Introduction 

The first deep survey above 20 keV performed by the 
INTEGRAL satellite has allowed the detection of more than 400 



* Based on observations obtained with XMM-Newton and 
INTEGRAL, ESA science missions with instruments and contri- 
butions directly funded by ESA Member States and NASA, and with 
Suzaku, a Japan's mission developed at the Institute of Space and 
Astronautical Science of Japan Aerospace Exploration Agency in 
collaboration with U.S. (NASA/GSFC, MIT) and Japanese institutions. 



X-ray sources dBirdetalJl2007l) . The extensive survey of the 
Galactic plane has revealed that the contribution of Galactic X- 
ray binaries, especially the c ataclysmic variable ( CV) type, is 
non-negligible in hard X-rays (ISazonov et al.ll2006l) . Most of the 
hard X-ray emitti ng CVs were found to be magnetic intermedi- 
ate polars (IPs) (Barlo w et alj|20 0^. Up to now, IPs represent 
~ 5% of the INTEGRAL sources detected at that energy range 
and this number is prone to increase in the near future, as demon- 
strated by systematic optical follow-ups f or some of the ~ 10 
INTEGRAL unidentified sources (see e.g. iMasetti et alj|2008h . 



2 



G. Anzolin et al.: Broad-band properties of the hard X-ray CVs IGR J00234+6141 and IRXS J213344.1+510725 



From the detection of thousands of discrete low-luminosity X- 
ray sources in a deep Chandra survey of regions close to the 
Galactic center, Muno et al. (2004) proposed that magnetic CVs 
of the IP type may represent a significant fraction of the Galactic 
background. Further support that magnetic CVs could still be a 
hidden population of faint X-ray sources and, therefore, play an 
important role in the X-ray emission of the Galaxy comes from 
recent studies of t he Galactic ridge with INTEGRAL and RXTE 
dRevnivtsev et al.ll2008l) . 

IPs are believed to harbor weakly magnetized (B < 10 MG) 
white dwarfs (WDs) because of their fast asynchronous rotation 
with respect to the orbital period and of the lack of significant 
polarized emission in the optical/near-IR in most systems. This 
contrasts with the other group of magnetic CVs, the polars, that 
instead possess strongly magnetized (B ~ 10 - 230 MG) WDs 
rotating synchronously with the binary period. The X-ray prop- 
erties of magnetic CVs are strictly related to the accretion mech- 
anism onto the WD primary. Material from the late type compan- 
ion is driven by the magnetic fie ld lines onto the magnetic polar 
caps, where a shock develops dAizul [1973 1) below which hard 
X-rays and cyclotron radiation are emitted. Bremsstrahlung ra- 
diation is believ ed to be the dominant cooling mechanism in IPs 
jWu et alJI 19941) . while cyclotron radiation may dominate in po- 
lars. The complex interplay between the two mechanisms greatly 
depends on both magnetic field strength and local mass accretion 
rates ( Woelk & Beuermann 1996; Fischer & Beuermann 2001). 
Hence, if IPs indeed host weakly magnetized WDs with respect 
to polars, this could qualitatively explain why they are hard X- 
ray sources. The detection of a soft X-ray optically thick com- 
ponent in an increasing number of systems jAnzolin et alj|2008l) 
poses further questions in the interpretation of the X-ray emis- 
sion properties of IPs. 

In the framework of an ongoing optical identifica- 
tion program, we identified two new members of the IP 
group: IRXS J213344. 1-H510725 ^ IGR J21335-H5105 
(hereafter RXJ2133) jBonnet-Bidaud et alj |2006|) and 
IGR J002 34-H614I = IRXS J0022 58.3-I-6141I1 (hereafter 
IGR0023) (lBonnet-Bida udetalJl20 07h Both of theni are hard 
CVs in the INTEGRAL source catalog dBird et alj|2007h . 

The weak hard X-ray source IGR0023 was detected by the 
INTEGRAL satellite during an observation o f the Cassiopeia re- 
gion of the Galaxy jden Hartog et al. 20061). The o ptical coun- 
terpart of IGR0023 was identified bv lMasetti et alj ( f2006) . who 
proposed a possible magnetic nature of this CV. A tentative 
~ 570 s optical pe riodicity was recognized in the R band 
jBikmaev et alj|2006l) . However, a clear periodic modulation of 
56 3.53 + 0.62 s was di s covere d with optical photometric data 
bv lBonne t-Bidau d et al.l (l2007l) and readily ascribed to the rota- 
tional period of the WD, while an orbital period of 4.033 + 0.005 
hr was derived from optical spectroscopy. The properties of the 
INTEGRAL spectrum, which is well fitted by a bremsstrahlung 
with a temperature of 3 1 ke V, strongly support the magnetic na- 
ture of IGR002 3. 

RXJ2133 was identified as a h ard X-ray point sou rce from 
the ROSAT Galactic Plane Survey dMotch et al.lll998h . A clear 
persistent optical light pulsation at 570.823 + 0.013 s was then 
discovered with fast optical photometry, while optical spec- 
troscopy revealed an additiona l periodic variability at 7.193 + 
0.016 hr jBonnet-Bidaud et al.l 12006) . These two periodicities 
were respectively identified as the WD spin and the orbital peri- 
ods, thus suggesting RXJ2133 as a member of the IP class with a 
relatively long orbital pe riod which falls into the so-called IP gap 
betwee n 6.5 and 9.5 hrs dSchenker et al.ll2004l) . lKataiainen et al.l 
(l2007h found also that RXJ2133 emits optical circularly polar- 



ized light up to ~ 3% and proposed that the WD magnetic field 
could be as high as 25 MG, one of the highest amongst IPs. 

The X-ray variability and broad-band spectra of these two 
systems have been investigated using pointed XMM-Newton 
jjansen et al.ll200l]) observations and publicly ava ilable hard X- 
ray d ata obtained with the INTEGRAL satellite jWinkler et all 
2003). In the case of RXJ2133, we also present the temporal and 
spectral analysis of a pointed Suzaku (.Mitsuda et al.,.2007.) ob- 
servation. 



2. Observations and data reduction 

The summary of the observations of IGR0023 and RXJ2133 is 
reported in Table [T] 



2.1. The XMM-Newton observations 

For all our XMM-Newton observations, we reprocessed a nd an- 
alyzed the E PIC-pn ( Struder et al. 2001), MOS ( T urner et al 
2001), RGS dden Herder etal.1 [2001]) and OM dMason et al 
2001) data using the standard reduction pipelines included in 
SAS 8.0 and the latest calibration files. For IGR0023, because 
of a problem with the standard source detection task, the OM- 
UVM2 data were reprocessed at MSSL using an unreleased 
reduction routine. Heliocentric corrections were applied to the 
EPIC and OM data of both sources. The SAS tasks rmfgen and 
arfgen were used to generate the photon redistribution matrix 
and the ancillary region files for all the EPIC cameras and RGS 
instruments. 

IGR0023 was observed on July 10, 2007 (OBSID: 
0501230201) with the EPIC-pn and MOS cameras operated in 
full frame imaging mode with the thin and medium filters, re- 
spectively. The total exposure times were 25 ks for EPIC-pn and 
26.6 ks for both the MOS cameras. The RGS was operated in 
spectroscopy mode for a total exposure time of 26.9 ks. The 
OM was operated in fast imaging mode using sequentially the 
V (5000-6000 A) and UVM2 (2000-2800 A) filters, for 9.8 ks 
each. 

A 28" aperture radius was used to extract EPIC light curves 
and spectra from a circular region centered on the source and 
from a background region located on the same CCD chip where 
the source was imaged. In order to improve the S/N ratio, we 
filtered the data by selecting pattern pixel events up to double 
with zero quality flag for the EPIC-pn data, and up to quadruple 
pixel events for the EPIC-MOS data. The average background 
level of the EPIC cameras was quite low for almost all the dura- 
tion of the observation, with the exception of a moderate flaring 
activity, which occurred during the last ~ 2000 s of the EPIC-pn 
exposure. This flare did not significantly affect the data used for 
the timing analysis, however we conservatively did not consider 
these events in the spectral analysis. 

Due to the weakness of IGR0023, the RGS data had poor 
S/N ratios and therefore were not useful for spectral analysis. 

Background subtracted OM-V light curves were obtained 
with a binning time of 20 s, while the OM-UV Ml light curves 
were provided by MSSL with a binning time of 10 s. The 
average count rates were 3.00 counts s"' in the V band and 
0.18 counts s~' in the UVM2 band, corresponding to instrumen- 
tal magnitudes V - 16.8 and UVM2 - 17.6 and average fluxes 

of7.4xl0""'erg cm"^ s"' A"' and4.0xl0""'erg cm"^ s"' A"', 
respectively. As a comparis on, the continuum flux of the op- 
tical spectrum obtained by iMasetti et alj (l2006i) was ~ 6 x 



10 erg cm ^ s 'A 
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Table 1. Summary of the X-ray observations of IGR0023 and RXJ2133. We report total exposure times for the INTEGRAL obser- 
vations and net exposure times for the others. 



Object 


Instrument 


Date 


UT (start) 


Exposure time (s) 


Net count rate (counts s ') 


IGR0023 


EPIC-pn 


2007-07-10 


05:58 


24660 


0.984 + 0.007 




EPIC-MOS 




05:36 


26540 


0.345 ± 0.004 




RGS 




05:35 


26706 


0.038 + 0.001 




OM-V 




05:44 


1960 


3.26 + 0.04 








Uo. ZZ 


1 O'^O 

i yj\) 


D.VZ ± U.U4 








07:01 


1960 


2.79 ± 0.04 








08:12 


1962 


2.97 + 0.04 








08:50 


1959 


2.85 + 0.04 




0M-UVM2 




09:28 


1960 


0.16 ±0.01 








10:06 


1959 


0.19 + 0.01 








10:45 


1960 


0.20 + 0.01 








11:23 


1960 


0.18 + 0.01 








12:01 


1960 


0.18 ±0.01 




IBIS/ISGRI 






~ 6900000 


0.10 + 0.01 


RXJ2133 


EPIC-pn 


2005-05-29 


13:33 


13 585 


5.42 ± 0.02 




EPIC-MOS 




12:47 


16 580 


1.398 ± 0.008 




RGS 




13:28 


13 890 


0.147 + 0.004 




OM-B 




10:11 


2 320 


24.9 + 0.1 








iU. JJ 


z jiy 


zj.z + U. 1 








11:39 


2319 


22.4 + 0.1 








12:23 


2319 


22.1+0.1 








13:07 


2 320 


23.7 + 0.1 




0M-UVM2 




13:51 


2319 


0.69 + 0.02 








14:35 


2 321 


0.68 + 0.02 








15:20 


2 320 


0.89 + 0.02 








16:03 


2 320 


0.85 + 0.02 








16:48 


2319 


0.87 + 0.02 




EPIC-pn 


2005-07-06 


18:03 


9 871 


5.12 + 0.03 










Vj OoU 






RGS 




17:04 


13910 


0.120 + 0.005 




0M-UVM2 




17:31 


1681 


0.60 + 0.02 








18:16 


1679 


0.70 + 0.02 








19:20 


1680 


0.64 + 0.02 








19:53 


1680 


0.69 + 0.02 








20:27 


1679 


0.59 + 0.02 




IBIS/ISGRI 






~ 3 740 000 


0.55 + 0.02 




XIS-FI 


2006-04-29 


06:50 


84288 


0.738 + 0.002 




XIS-BI 




06:50 


84288 


0.892 + 0.004 




HXD 




06:50 


62 879 


0.739 + 0.004 



RXJ2133 was observed on May 29, 2005 (OBSID: 
0302100101) with all the EPIC cameras operated in full frame 
mode and with the medium filter Due to high background ra- 
diation the observation time was shortened and the net exposure 
times were 13.6 ks for EPIC-pn, 16.5 ks for EPIC-MOS and 13.9 
ks for RGS, the latter operated in spectroscopy mode. The OM, 
operated in fast imaging mode, was used sequentially with the B 
filter, covering the spectral range 3900^900 A, and the UVM2 
filter for 1 1.6 ks each. RXJ2133 was observed again on July 06, 
2005 (OBSID: 0302100301) with the same instrumental config- 
urations and with net exposure times of 9.9 ks for EPIC-pn, 13.7 
ks for EPIC-MOS and 13.9 ks for RGS. The OM was only op- 
erated with the UVM2 filter for a total exposure time of 8.4 ks. 
The source was found at about the same count rate at the two 
epochs in all instruments. 

EPIC light curves and spectra of source and background 
were extracted from circular regions of 37" radius. The same fil- 
tering adopted for IGR0023 was applied to improve the S/N ratio 



for both cameras. During the observation of May 2005 the back- 
ground was moderately active, but not enough to affect the tim- 
ing analysis. Filtering of higher background periods was done 
only during the extraction of spectra from both EPIC and RGS 
data. These spectra were rebinned to have a minimum of 25 and 
20 counts per bin, respectively. 

OM B and UVM2 light curves were extracted with a bin- 
ning time of 10 s and 20 s, respectively. The average count 
rates were about 25 counts s"' in the B band and 0.7 counts s"' 
in the UVM2 band, corresponding to instrumental magni- 
tudes B = 15.8 and UVM2 = 16.2. These translate into B 
and UVM2 fluxes of 3 x 10"'^ erg cm"^ s"' A and 1.5 x 
10"'^ erg cm"^ s"' A , respectively. 

2.2. The INTEGRAL observations 

The INTEGRAL IBIS/ISGRI dUbertini et al.ll2003l: iLebrun et al.1 
2003i) hard X-ray data of both sources were extracted from all 



4 



G. Anzolin et al.: Broad-band properties of the hard X-ray CVs IGR J00234+6141 and IRXS J213344.1+510725 



pointings within 12° from the source positions, spanning from 
March 2003 to October 2006. The total effective exposure times 
are ~ 6.9 Ms (2875 pointings) and ~ 3.74 Ms (1563 pointings) 
for IGR0023 and RXJ2133, respectively. To study the weak per- 
sistent X-ray emission, the time averaged ISGRl spectra have 
been obtained from mosaic images in five energy bands, loga- 
rithmic spaced between 20 and 100 keV. Data were reduced with 
the standard OSA software version 7.0 and, then, an alyzed using 
the algorithms described bv lGoldwurm et akl ( l2003l) . 

2.3. The Suzaku observations ofRXJ2133 

RXJ2133 was observed with Suzaku between Apr 29, 2006 
and May 1, 2006 (sequence number 401038010). We have 
analyzed data from the X-ray imaging spectrometer (XIS, 
iKovama et al. J 20071)) and the no n-imaging hard X-ray detec- 
tor rHXD. iTakahashi et al.1 (I2007h 1. The observation was done 
with the object at the "HXD-nominal" pointing position, ~ 5' 
off-axis from the center of field-of-view (FOV) of the XIS, to 
optimize the S/N ratio of the HXD data. We based our analysis 
on data processed using the V2.0.6.13 pipeline released as a part 
of HEADAS 6.3.1. 

For the XIS, we updated the energy scale calibration using 
the February 1, 2008 release of the calibration database. We then 
applied the following screening criteria: attitude control system 
in nominal mode, pointing within 1'5 of the mean direction, 
XIS data rate medium or high, the satellite outside the South 
Atlantic Anomaly (SAA) and at least 180 s after the last SAA 
passage, elevation above Earth limb > 5°, elevation about the 
bright Earth limb > 15°. An inspection of the XIS image re- 
vealed a second source, near the center of the FOV, with a flux 
of 6.7 X 10"^^ ergcm"^ s ' in the 2-10 keV band. This same 
source is seen in the XMM-Newton data at similar flux levels Q. 
Although being faint, we conservatively excluded this source 
from our analysis of the XIS data. We cannot do so for the HXD 
data analysis, but it is expected to have a negligible effect. 

We used a 3 '5 radius extraction region and an annular extrac- 
tion region of 7'5 outer and 4' inner radii for the background, 
both centered on the position of RXJ2133. For spectroscopy, we 
summed the data and the responses of the three XlS-Fl units be- 
cause they have nearly identical responses. For photometry, we 
added background subtracted light curves from all the XIS units 
over the energy range 0.3 (Fl) / 0.2 (Bl) - 12 keV. 

For the HXD data, we took the PIN event data from the pro- 
cessing pipeline and applied the dead time cor rection. We ob- 
tained the "tuned" non X-ray background files (Fukaza wa et alj 
|2009|), estimated by the HXD team using LCFITDT method. 
For phase-resolved spectroscopy, we used the phase-averaged 
background and dead-time fraction, since both tend to vary on 
a longer time scale. 



3. Data analysis and results 

3.1. IGR0023 

3.1 .1 . The X-ray variability 

The EPlC-pn light curve extracted in the energy range 0.2-10.0 
keV and binned in 20 s time intervals reveals a clear variabil- 
ity with a time scale of the order of ~ 10 min. We have not 
analyzed the EPIC-MOS light curves because they were too 
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' We have not investigated the nature of this second source, but a 
likely explanation is that it is a background AGN. 



Fig. 1. Power spectra of the EPlC-pn data of 1GR0023 in se- 
lected energy ranges. From bottom to top: 0.2-10.0 keV, 0.2-1.0 
keV, 1.0-2.0 keV and 2.0-10.0 keV. Vertical dotted line indicate 
the frequencies w, 2a) and fi (see text). 



noisy to provide reliable results. The power spectrum of the full- 
band EPlC-pn light curve (see Fig. [1) shows significant peaks 
at (L) ~ 154 d"', at 2a) and also at low frequencies. A peak at 
fl ~ 22 d"', although its significance is b elow 2cr, is close to a 
pseudo-periodicity detected in the optical (iBonnet-Bidaud et al.l 
i2007h . 

A sinusoidal fit was performed on the EPlC-pn light curve 
previously corrected for low frequency trends by using a third- 
order polynomial. We used three sinusoids with different fre- 
quencies accounting for all the observed peaks, thus finding 
w= 153.83+0.15d"' and/i = 21.94+0.15 d"' (errors are at the 
Icr confidence level). The inferred period P^j = 561.64 + 0.56 s 
can be identified with th e WD spin period, since the difference 
with the optical period of lBonnet-Bidaud et al. I (l2007h IS not sig- 
nificant. 

A Fourier analysis was also performed on EPlC-pn light 
curves extracted with a 40 s binning time in selected energy 
bands: 0.2-1.0 keV, 1.0-2.0 keV and 2.0-10.0 keV (see Fig.fl. 
Peaks at a) and 2a> are clearly detected in the low and interme- 
diate energy bands, while they do not appear at high energies. 
Instead, the pseudo-periodicity at low frequency is visible only 
in the intermediate 1-2 keV band. 

We then folded the EPIC-pn light curves with the 561.64 
s spin period (Fig. |2]i using the time of maximum obtained 
from the sinusoidal fit of the 0.2-10.0 keV light curve: HJD = 
2454291.8667(2). These light curves show a quasi-sinusoidal 
modulation only below 2.0 keV, with a secondary maximum at 
<p ~ 0.5. The pulse amplitude is ~ 50% in the 0.2-1.0 keV range 
and ~ 16% in the 1.0-2.0 keV range. The count ratio between 
the 1.0-2.0 keV and the 0.2-1.0 keV bands indicates an harden- 
ing of the emission at the secondary maximum. 
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Fig. 2. The spin pulse of IGR0023 in selected energy ranges of 
the EPIC-pn and in the OM-V band. All the light curves are 
folded at the 561.6 s X-ray spin period using the time of max- 
imum quoted in the text. Solid lines represent sinusoidal fits to 
the corresponding data. 



Fig. 3. The EPIC-pn and INTEGRAL spectra of IGR0023 are 
shown with model B discussed in the text and parameters re- 
ported in Table |4] (we do not report the combined MOS spec- 
trum for clarity). Each spectral component that contributes to 
the best-fit model is shown separately. The bottom panel shows 
the residuals expressed in terms of cr. 



3.1 .2. The visible and UV liglit curves 

The power spectrum of the OM-V light curve shows a strong 
peak at the X-ray spin period, while that of the UVM2 light 
curve does not reveal any significant peak. The average count 
rate is not constant in the 5 observations with the V filter span- 
ning ~ 3/4 of the orbital cycle, probably suggesting a depen- 
dence on the orbital period. The spin-folded OM-V light curve 
(lower panel of Fig|2]) is pretty sinusoidal, with an amplitude 
of 8 + 1 %, and shows a single peak at a phase consistent with 
that of the main maximum of the X-ray pulse. We notice that 
iBonnet-Bidaud et"an ( |2007|) also found a single peaked pulsa- 
tion in their observations carried out with the Gunn-g filter 

3.1.3. Spectral properties of IGR0023 

The EPIC-pn and combined MOS spectra (in the range 0.3-10.0 
keV) and the IBIS/ISGRI spectrum (in the range 20-100 keV) 
were simultaneously analyzed with the XSPEC 12 package. An 
absorbed isothermal optically thin MEKAL component plus a 
zero width Gaussian at 6.4 keV fit relatively well the spectrum 
(xl - 1.04), but the temperature is unconstrained (> 78 keV). 
The fit improves using a multi-temperature CEMEKL emission 
component and a dense {Nh ~ lO^"* cm"^) absorber covering 
~ 40% of the source (model A in Table |4]i, although the a pa- 
rameter assumes an unreasonably high value. Metal abundances 
are consistent, within errors, with the solar values. The total ab- 
sorber is likely of interstellar origin, as it is comparable to that 
of the total Galactic ab sorption in the direction o f the source 
(A^H.gai = 7.4x 10^' cm^. lDickev & Lockmanl ( ll990l) ). The dense 
partial absorber, instead, is likely located close to the source, as 
suggested by the energy dependence of the spin light curve. We 
also obtained similar quality fits by substituting the CEMEKL 



component with 2 (Fig. O or 3 MEKALs (model B and C in 
Table m respectively). In both cases, we find a low (kT ~ 0.17 
keV) and an intermediate temperature (kT ~ 10 keV) compo- 
nents, and, in model C, a lower limit to the temperature of the 
third MEKAL kT > 27 ke V. The observed X-ray flux in the 0.2- 
10.0 keV range is 6.8 x 10"'^ erg cm"^ s"'. We notice that only 
model B and model C account for the O VII (21.9 A) and O VIII 
(19.1 A) lines barely detectable in the low quality RGS spectra. 

The INTEGRAL spectrum above 60 keV seems to be under- 
predicted by the three best-fit models but, given the poor S/N 
ratio, the excess of counts is not significant. The addition of 
a reflection component, suggested by the presence of the Fe 
line at 6.4 keV, does not improve the fits. We also fitted the 
broad-band spectru m of IGR0023 with the rec ent post-shock re- 
gion (PSR) model of Suleim anov et all (12008') that computes the 
emergent spectrum taking also into account the Compton scat- 
tering (x^ I d.o.f. = 136/110). We obtained a shock temperature 
^rshock = 51 + 11 keV and absorption parameters consistent with 
those found with models A, B and C. The contribution of the 
Compton scattering is found to be ~ 10% and, within errors, 
does not seem to affect the temperature determination. 

A phase-resolved analysis of the EPIC-pn spectrum of 
IGR0023 was performed selecting two phase intervals centered 
on the two maxima of the spin pulse (0 - 0.85 - 1.15 and 
4> - 0.35 - 0.65, respectively). We used model B, since model A 
and C would give badly constrained parameters. The hydrogen 
column density of the total absorber, the temperatures of the two 
MEKALs and the metal abundance were kept fixed to the val- 
ues obtained for the average spectrum. As shown in Tabled the 
normalization of the low-temperature optically thin component 
is significantly lower at the secondary maximum, where we also 
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Table 2. Spectral parameters obtained from fitting the EPIC-pn 
spectra of IGR0023 extracted in the phase intervals quoted in the 
text with model B shown in Table |4] 



Parameters 


Maximum 1 


Maximum 2 


A'h (10"' cm-') 


1 5+'" 

-0.6 




Gov. Frac. 


"•-'"-0.07 


0.51 ±0.04 


C, (10-^) 


1 2+"-3 


0.5 ±0.5 


C2 (10-2) 


^•-^-0.1 


1.7 ±0.1 


-f"o.2-10.0 kcV 


2.06 


2.00 


(10-'2 erg cm-- s-') 








1.01 (259/257) 1.10(258/235) 
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Fig. 4. Power spectra of the Suzaku XIS light curves of RXJ2133 
in selected energy ranges. From bottom to top: 0.2-12.0 keV, 
0.2-0.8 keV, 0.8-1.2 keV, 2.0-4.0 keV and 6.0-12.0 keV. 
Vertical dotted lines indicate the frequencies oj and lu (see text). 
The spectral window is shown in the inset. 



find a marginal evidence of an increase of the covering fraction 
of the local absorber. 



3.2. RXJ2133 

3.2.1. The X-ray periodicities in RXJ2133 

The EPIC light curves, extracted in the 0.2-12.0 keV energy 
range and binned in 5 s time intervals, and the Suzaku XIS light 
curve, extracted with a binning time of 8 s, show a variability of 
the order of 10 min. In the EPIC light curve of May 2005 a quasi- 
sinusoidal trend with an apparent period of ~ 2 hr is found, while 
in July 2005 a ~ 1 hr pseudo-periodicity is detected. However, 
in the later XIS observation we do not find any evidence of such 
a variability. 

We then performed a Fourier analysis on the light curves at 
the three epochs. Peaks at the optically identified spin frequency, 
ll), and at 2a) are clearly detected. In Fig. |4]we report the full- 
band XIS power spectra, as well as that obtained in different 
energy bands. Sinusoidal fits to the full-band light curves give 
OJ = 151.30+0.25d"', w = 151.55+0.78d"' andw = 151.350 + 
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Fig. 5. Spin-folded light curves of RXJ2133. From top to bottom: 
full-band XIS of April 2006, EPIC-pn and 0M-UVM2 of July 
2005, EPIC-pn and OM-B of May 2005. Sofid lines represent 
sinusoidal fits to the corresponding data. 



0.009 d-' for the May and July 2005 XMM-Newton and April 
2006 Suzaku data sets, respectively. In the analysis of the EPIC- 
pn light curves, a third sinusoid has to be included, accounting 
for the low frequency variations, thus giving /May2005 = 43.78 + 
0.29 d"' and /juiy2005 = 27.87 + 0.41 d"' . We identify the precise 
Suzaku period = 570.862 + 0.034 s with the true spin period 
of the accreting WD. This agrees, within errors, with the optical 
570.823 + 0.013 s period. 

We folded the EPIC-pn and XIS fight curves at the 570.86 
s X-ray spin period using the time of maxima of the pulsa- 
tion at (jj: //yZ)May2005 = 2453520.1444(2), HJDMyims = 
2453558.3134(4) and HJDApdmm = 2453856.09618(5). The 
three folded light curves (Fig.|5]l present two maxima at phases 
0.9 and 0.35, with a dip at phase ~ 0.1 that is more evident in 
the May 2005 data. The full amplitude of the primary maximum 
at <p ~ 0.9 is almost similar at the three epochs (~ 36%), while 
that of the secondary maximum decreases from ~ 13% in May 
2005 to ~ 3% in 2006. 

The energy -resolved EPIC-pn light curves (Fig.|6]l generally 
show the secondary maximum dominating the emission in the 
soft 0.2-0.5 keV band and the primary maximum dominating in 
the range 0.5-2.0 keV. At higher energies the two maxima have 
similar amplitudes of ~ 10%. The XIS light curves, though being 
broadly similar to the EPIC-pn light curves, present an evident 
secondary maximum only in the hard X-rays. 

3.2.2. line UV and optical liglit variability of RXJ2133 

The power spectrum of the OM-B light curve shows a strong 
peak at the spin frequency (consistent within errors with the 
X-ray and the previous optical determinations). We also detect 
power at low frequencies, but we cannot establish whether this 
indicates a true periodic variability. The OM-B light curve of 
May 2005 folded at the spin period (see Fig.|5]l is single peaked. 
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Fig. 6. EPIC-pn (May 2005 and July 2005) and XIS spin-folded 
light curves of RXJ2133 extracted in selected energy ranges. All 
the light curves are folded using the time of maximum quoted in 
the text. 



with a broad maximum at phase ~ and a full amplitude of 
14%. The UV light curve, instead, is almost unmodulated in May 
2005, while in July 2005 it is sinusoidal with a maximum shifted 
by half a period with respect to the OM-B light curve of May 
2005. A similar anti phased behavior is also seen in some soft X- 
ray IPs, like UU Col (Ide Martino et al.l2006l) . However, we point 
out that each OM light curve covers only ~ 41% of the orbital 
period and hence could be affected by orbital dependent changes 
of the spin modulation. 

3.2.3. The X-ray spectrum of RXJ2133 

The X-ray spectrum of RXJ2133 was studied in the broad-band 
energy range 0.3-100 keV using the EPIC-pn and combined 
MOS data at each epoch, together with the IBIS/ISGRI spec- 
trum. We also analyzed the combined XIS -FI, XIS-BI and HXD 
spectra in the restricted range 0.3-40 keV, where the S/N ratio is 
still good. The average X-ray spectrum is well fitted by a multi- 
temperature optically thin component plus a rather hot {kTBs ~ 
100 eV) blackbody, that is required to best fit the data at ener- 
gies below 1 keV. A total absorber with a low hydrogen column 
density, likely of interstellar origin (the total Galactic absorption 
along the direction of the source is A^H.gai = 7.35 x 10^' cm"^, 
iDickev & Lockmanid 19901) ). and a dense (Nh ~ 10^^ cm ^) ab- 
sorber partially covering the source are also required. In all the 
spectral fits we fixed the Nh of the total absorber to the value 
found for the May 2005 data, since it is not expected to vary 
with time, and used an unresolved Gaussian line to account for 
the fluorescence Fe line at 6.4 keV. The best-fit parameters are 
reported in Table|5] where the optically thin components in mod- 
els A, B and C are the same used in the corresponding spectral 
models of IGR0023. We have marked with 1 and 2 the May and 
July 2005 XMM-Newton spectra, respectively, while 3 indicates 



IRXS J2133. 7 + 5107 

-I I ^ 

EPIC-pn 



Folded model 

ckbody 
MEKAL (1.1 keV) 
MEKAL (4.B keV) - 
MEKAL (39 keV) 
Gaussian 



IBIS/ISGRI 




10 

Energy (keV) 
Uglll WeXHgK FeXVII 



l.OE-02 




lO.O 



15.0 30,0 
tfavdeiigth [S] 



35,0 



Fig. 7. Upper figure: the EPIC-pn (May 2005) and INTEGRAL 
spectra of RXJ2133 are shown with model C discussed in the 
text and parameters reported in Table |5] (the combined EPIC- 
MOS spectrum is not plotted for clarity). The different spectral 
components contributing to the best-fit model are shown sepa- 
rately. Bottom panel shows the residuals expressed in terms of 
cr. Lower figure: the RGS spectra of RXJ2133. 



the Suzaku spectra 0. Also for this object, the INTEGRAL spec- 
trum is underpredicted above ~ 60 keV. 

In general, we find no significant differences between the 
spectral parameters of May 2005 and July 2005, with the ex- 
ception of the normalization of the blackbody component that 
is slightly higher in the earliest epoch. The absorbed flux in the 
range 0.2-10.0 keV is 2.4 x 10"" ergcm"^ s"'. In the Suzaku 
spectra, the temperature of the hottest MEKAL (or the maxi- 
mum temperature of the CEMEKL model) attains higher val- 
ues than in 2005. In addition, we find a decrease of the cov- 
ering fraction of the partial absorber (50%) and a lower value 
of the normalization of the blackbody (with the exception of 



- We caution that the quoted errors reflect only the statistical errors. 
However, our experimentation of adjusting the HXD non X-ray back- 
ground by ± 1 % from the best estimate values did not cause significant 
changes in the fit parameters, showing that the systematic errors are not 
dominant in this case. 
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Fig. 8. The Suzaku spectra of RXJ2133 are shown with model 
C discussed in the text and parameters reported in Table |5] The 
bottom panel shows the residuals. 



model C). The absorbed flux in the 0.2-10.0 keV range is 
2.5 X 10"" erg cm"^ s"', in agreement with that found in the 
two observations of 2005. Although we find diff'erent values of 
the bolometric fluxes of both the blackbody and the optically 
thin components between the two epochs, they are consistent 
within errors: Fbb = (4.1:^o4 - 4.8+J}g) x 10"" erg cm"^ s"' and 

r+0.8N ^ in-ii 



i^hard = (7.3!|:i - T.e^O;^) X 10-" erg cm-^ s"' in 2005, Fbb 
4.0+^ ° X 10"" erg cm^^ s"' and f hard > 8.3 x 10"" erg cm^^ s"' 
in 2006. Therefore, the ratio between the soft and hard X-ray 
fluxes is of the order of 0.5. 

Because of the similar temperatures of the optically thin 
components of model C found at the three epochs, we adopt 
this model to describe the broad-band spectrum of RXJ2133 (see 
Fig.Qand Fig.[8j. Although there is clear evidence of a temper- 
ature gradient in the post-shock region, the broad-band spectra 
do not provide constraints on lower temperatures. Instead, the 
presence of a low temperature plasma can be inferred from the 
May 2005 RGS spectra, where the strong O VIII (0.651 ± 0.002 
keV , E.W. = 22 eV) and NeX (1.02 + 0.01 keV, E.W. = 14 
eV) lines are clearly visible (see Fig. |7]i. The O VII line at 0.58 
keV {E.W. = 8 eV) and the Fe XVII line at 0.73 keV {E.W. = 1 1 
eV) are likely present, while there is only a weak evidence of 
the Ne IX line. We then estimate a ratio ~ 2.5 between the H- 
and He-like oxygen lines, that would imply a temperature of 
~ 0.3 - 0.4 keV. This is also supported by the Ne line ratio of 
~ 1.5, again indicating kT ~ 0.4 keV. 

The Fe line at 6.4 keV has a large equivalent width {E.W. - 
150 - 170 eV), suggesting a Compton reflection component. 
However, we do not find improvements of the fits with the inclu- 
sion of a reflection component. We then applied the PSR model 
to the combined XMM-Newton (May 2005) and INTEGRAL 
spectra {x^ / d.o.f. = 347/279), as well as to the Suzaku spec- 
tra {x^ I d.o.f. - 395/359). We find a shock temperature of 50 + 2 
keV and 53 + 3 keV, respectively, while the parameters of the ab- 



sorbers and of the blackbody component are consistent, within 
errors, with those found using models A, B and C. 

We also performed a phase-resolved analysis of the EPIC- 
pn (May 2005) and XIS spectra. The former was extracted 
in phase intervals approximately centered on the two maxima 
(0 = 0.72 - 0.94 and 4> = 0.16 - 0.28) and at the minimum 
{^ - 0.44 - 0.61). Instead, the higher quality XIS spectrum was 
extracted in 4 consecutive phase intervals, namely (p - 0.05-0.3 
(PI), 0.3-0.55 (P2), 0.55-0.8 (P3) and 0.8 - 1 .05 (P4). We fitted 
these spectra adopting model C, fixing Nh of the total absorber, 
the temperatures of all the emission components and Az to the 
values found for the corresponding average spectrum. The re- 
sults reported in Table |3] show significant differences in some of 
the parameters. The normalization of the soft MEKAL increases 
at the primary maximum in the two epochs, whilst only upper 
limits are found at other phases. Also, a decrease of the hydro- 
gen column density of the partial absorber and a slight increase 
of the blackbody normalization at the secondary maximum are 
found only in May 2005. 



4. Discussion 

4.1. IGR0023 

The XMM-Newton data confirm the previously identified opti- 
cal 561 s period as the spin period of the accreting WD. A 
pseudo-periodicity of the order of 1 hr, though with a signif- 
icance level below 3cr, seems to be present a bove 1 keV and 
was previously found also in optical data (Bon net-Bidaud et~aLl 
'2007). This variability, if real, is not easy to explain. Assuming a 
circular Keplerian orbit, a periodicity of 1 .09 hr implies a radius 
of ~ 3 x lO'" cm. For comparison, adopting a mass ratio q - 0.5 
and a WD mass of ~ 0.9Mq (see below), the outer edge of the 
accretion disc would be 4 x 10'" cm. This may suggest an origin 
within the disc, but unlikely due to ab sorption effects as instead 
found in many IPs dParker et al.ll2005h . 

The spin pulse is detected only below 2 keV, where the 
power spectrum shows substantial signal also at 2a). The double- 
humped spin pulse and the increase of amplitudes with decreas- 
ing energy could be the signatures of two accreting poles, as well 
as of variable absorption. From the phase resolved spectra, we 
find a decrease by a factor ~ 2 in the normalization constant of 
the MEKAL at 0. 17 keV, as well as an increase of a dense partial 
absorber at the secondary maximum. This prevents us to isolate 
the true contribution from a secondary pole. In addition, the opti- 
cal pulse is single-humped, thus not revealing a secondary pole. 
Hence, we are left with the uncertainty on whether the X-ray 
pulse shape is due to the presence of a secondary and less active 
pole or to phase-depende nt absorption effects in the accretion 



curtain jRosen et 



de pende i 
all 19881) 



The X-ray emission in IGR0023, extending up to ~ 90 keV, 
requires a multi-temperature optically thin plasma with a tem- 
perature distribution more complex than a simple power law. The 
spectral fits indicate a likely maximum temperature of ~ 50 keV. 
High plasma temperatures were also inferred from broad-band 
analysis of the combined XMM-New ton - INTEGRA L spectra 
of IRXS J173021.5-055933 (de M artino et al.ll2008l) and from 
hard X-ray observations of several other IPs ( Landi et ani2009l: 



ISuleimanov et al.l [20081; iBrunschweiger et aT 



20091). Alfliough 



caution has to be taken to interpret these temperatures in terms 
of the WD mass, the extension of spectra above 30 keV allow 
us to directly observe the exponential cut-off of the underlying 
continuum and hence the inferred temperatures are likely to be 
more reliable than those based only on softer energy coverages 
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Table 3. Spectral parameters obtained from fitting model C shown in Table |5] to the EPIC-pn (May 2005) and XIS spectra of 
RXJ2133 extracted in the phase intervals quoted in the text. 





EPlC-pn (May 2005) 


XIS (April 2006) 


Parameters 


Maximum 1 


Maximum 2 


Minimum 


PI 


P2 


P3 


P4 


Cbb (10-^) 


4.7 ± 0.5 


0.6 ± 0.2 

56+°"* 
6 4+' ° 


1 9+0.6 

^-■^-0.4 

55+" "^ 

"■-^-^-0.07 

4 6+"-^ 


1.1 ±0.2 
0.44 ± 0.02 
4.7 ± 0.2 


1.0 ±0.1 
0.53 ± 0.01 
4.9 ±0.2 


1.1 ±0.2 
0.46 ± 0.02 

4.2 ±0.2 


1.1 ±0.2 
0.38 ± 0.02 
4.4 ± 0.2 


C, (10-") 


14 ± 3 


< 2.5 


< 4.2 


< 0.8 


< 0.2 


< 1.5 


6.4 ±0.8 


C, (10-3) 
C3 (10-2) 


< 2.3 
1.6 ±0.2 


5 2+2-7 
1.2 ±0.2 


2 6+2-5 
1.4 ±0.2 


3 6+" ** 
1 48+°"'' 


3 9*°-^ 

15g°J.07 

^•-'-'-0 Ofi 


2.7 ±0.7 
1 47+0.09 


T 7+0.6 
•'■'-0.7 
1.32 ±0.08 


-f 0.2-10.0 kcV 


2.51 


2.38 


2.12 


2.52 


2.48 


2.35 


2.47 


(10-" ergcm-2 s"*) 
















xi 


1.04 


0.96 


0.97 


0.82 


0.87 


0.80 


0.87 


(x'/d.o.f.) 


(434/418) 


(266/276) 


(280/290) 


(1651/2019) 


(1721/1982) 


(1548/1939) 


(1747/2013) 



(e.g. XMM-Newton). For IGR0023 we est imate a WD mass of 
0.91+n |g/Vffri, consistent with that found bv lBrunschweiger et alj 
( |2009|) from the analysis of Swift/BAT data. Although relatively 
massive WD primaries were inferred in several ma gnetic CVs 
(see e.g. iRamsavl l200(il : iBrunschweiger et al.l I2009X the mass 
distribution of magnetic CV primari es is not very m uch different 
from that of non-magnetic systems ' Ramsavl (l2000h . Hence, for 
a CV to be a hard X-ray source, there should be an additional 
parameter that plays an important role. 

The hard X-ray bolometric flux Fhaid - 1.8 x 
10~" erg cm~^ s~' is an approximate estimate of the total accre- 
tion luminosity, as the pulsed components at other wavelengths 
need to be included as well as a secure distance estimate is miss- 
ing. Here we can obtain only a lower limit using the 8 % pulsed 
V band contribution, adopti ng E{B - V) - .35 as derived by 
the Nh of the total absorber (iRvter et al.ll 19751) . Also, assuming 
a lower limi t of 500 pc for th e distance of IGR0023, as esti- 
mated by Bo nnet-Bidaud et all (12007) using the 2MASS /Tband 
mag nitude when attributed solely to the secondary star (iKniggd 

I2OO6 ). we obtain Laccr > ^hai-d + Ly - 5.5 X 



10''4,0pcergs-'. 
This in turn gives us a lower limit for the accretion rate M > 



2.5 X 10 t^soopc g ^ ■ Unless the distance is much larger and 
the hard X-rays do not trace the bulk of accretion, this value is 
much l ower than 2.2 X 10'^ g s-' predicted for its 4.033 hr orbital 
period (IWarneiill995h . 

The short spin p eriods of IGR 0023 and RXJ2133 a re similar 
to tho se of YY Dra (iPatterson et al. 1992), V405 Aur dSkillmanl 
fT996h and IRXS J070407.9-H262501 ( Gansicke et al.l l2005h . 
Fast r otators were propo sed to possess weakly magnetized 
WDs jNorton et all 1 19991) . although V405 Aur and RXJ2133 
are among the few IPs showing polarized optical emissions. 
Assuming a pure disc accretion, as indicated by the absence of 
orbital sidebands in the power spectrum, a lower limit for the 
magnetic moment yU > 9x lO-'' c/soopc G cm"* is obtained using the 
inferred values of /Wwd and M. This could suggest that IGR0023 
is not strongly magnetized and not spinning at equilibrium. 



4.2. RXJ2133 

Our X-ray analysis confirms the 5 71 s period detected in op- 
tical photometry and polarimetry ( Bonnet-Bidaud et alj l2006t 
iKataiainen et al.ll2007l) as the rotational period of the WD pri- 
mary. RXJ2133 is, therefore, one of the most asynchronous sys- 
tems among IPs, with a ratio between the spin and the orbital 
periods of 0.022. 

The double peaked X-ray pulse shows a complex energy de- 
pendence and might be interpreted in terms of absorption and 



two emitting poles. The energy-resolved light curves suggest 
that the primary pole dominates the emission in the energy range 
0.5-5.0 keV, while the two poles contribute similarly at higher 
energies. This is also found in the phase-resolved spectral anal- 
ysis where the main changes are due to the emitting volume of 
an intermediate temperature region. The soft X-ray blackbody 
component is seen at all spin phases, suggesting either that the 
irradiated main pole is never occulted or that reprocessing oc- 
curs at both poles. However, in 2005 the secondary maximum 
is dominant in the soft band, and at this phase we find a de- 
crease in the absorption and a slight increase in the blackbody 
normalization. Hence, when the main pole points away from 
the observer, we are viewing the contribution of the irradiated 
area of the secondary pole. In order to have a lower absorption 
at this phase, a modificati on of the standard accretion curtain 
model dNorton et al.lll999h . where the optical depth is largest 
when viewing the curtain perpendicularly, could be envisaged. 

UV and optical data also support a scenario with two ac- 
creting poles whose contributions are however variable with 
time. The XMM-Newton UV observation of July 2005 shows 
a weak modulation antiphased with the X-ray pulse, while in 
May 2005 the UV light is not modulated but the B band pulse is 
single-peaked. Previous white light photometry acquired in 2003 
(Bonnet-Bidaud et al. 2006) shows the presence of two poles, 
being the first harmonic of both the spin frequency and the beat 
clearly detected. This is not seen in the 2006 optical photometry 
and polarimetry ( IKataiainen et alj|2007h as well as in the Suzaku 
X-ray data in the same year. Note that i Bonnet-Bidaud etaTI 
(2006) suggested a relatively low binary inclination (i < 45°) 
and the behavior of the circular polarization a long the spin cy- 
cle in dicates a magnetic colatitude /3 ~ 90° - ; (IKataiainen et al.l 
|2007|) . This implies ; > 45° and hence that the secondary pole, 
whenever active, can be observed. 

The broad-band X-ray spectrum of RXJ2133 is well fitted 
by a multi-temperature optically thin plasma with a likely max- 
imum temperature fcr^ax ~ 50 keV. As for IGR0023, a power- 
law multi-temperature flow appears a too simple and inadequate 
description of the post-shock emitting region, implying that the 
emergent spectrum is highly sensitive to local pressure and tem- 
perature across the flow. Furthermore, a blackbody component at 
100 eV is also required to account for the soft X-ray emission. To 
date this component has been detected in ~ 42% of the IPs, thus 
appearing a com mon characteristic s of magnetic CVs and not 
solely of polars (A nzolin et al.ll2008l) . If produced by the heating 
of a small area surrounding the WD polar cap, the range of val- 
ues of the normalization constant obtained by applying model C 
to the spectra of RXJ2133 (see Table|5]l implies an emission area 
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of 1.6 - 1.9 X lO'-' "^eoopc Here we use a minimum distance 
to the source as derived by iBonnet-Bidaud et"al] (|2006|) using 
2MASS K band magnitude. Assuming kT^^^^ ~ kT^bock, we in- 
fer a WD mass of 0.93 ± 0.04M(T), in agre ement with the values 
found by iBonnet-Bidaud et"aLl (l2006h and 'Brunschweig er et al] 
( 120091) . and a radius of 5.6 X 10^ cm ([Nauenberg 1972|). The 
fractional area of the blackbody emitting region is / ~ 2 - 
3 X 10"^ '^eoopc- ^^^^ could represent a small core onto the WD 
pole, as found in other soft X-ray IPs with similar temperatures 
(lHaberletaLll2002l: Ide Martino et alj|2008l) . We also point out 
that a ~ 100 eV blackbody would imply that the emission from 
the accretion region locally exceeds the Eddington limit. At that 
temperature, the radiation pressure is 4 times stronger than the 
gravity at the surface of a 0.93Mq WD. Clearly we still need a 
better understanding of emission mechanisms to solve this puz- 
zle. 

The bolometric flux of the blackbody represents ~ 40% 
of the total bolometric X-ray emission. RXJ2133 shows circu- 
larly polarized emission up to ~ 3% and hence the flux due 
to cyclotron radiation cannot be neglected. Since the blackbody 
flux does not exceed that of the hard X-ray emission, the re- 
processed radiation seen in the soft X-rays does not balance 
the primary radiation components. We have found that the UV 
emission is modulated by 6% at the spin period in July 2005, 
but antiphased with respect to the hard X-ray component, thus 
suggesting that the reprocessed radiation at the WD surface is 
also emitted at UV wavelengths. The pulsed UV flux observed 
in July 2005, dereddened for E{B - V) - 0.25, is Fuvm2 ~ 

2.1 X 10"'^ ergcm"^ s"' A . This value is much larger than 
that expected from the extrapolation of the soft X-ray flux to- 
wards UV wavelengths. Also, the UV flux integrated over the 
UVM2 band is only 4 % of the soft X-ray bolometric flux and 
hence only provide a lower limit to the UV luminosity. 

To evaluate the mass accretion rate we then consider the 
hard and soft X-ray components, as well as the modulated frac- 
tion of the flux in the B band of the May 2005 observation. 
A lower limit for the accretion luminosity is therefore: Laccr ^ 
5 X 10^3 erg s"'. We obtain M > 2.3 x 10'^ d^,^^ g s"', 

much lower than the secular value of 1.9 x 10'** g s"' predicted 
for its long 7.2 hr orbital period. From the detection of circu- 
lar polarization in the o ptical range, that peaks in the V band, 
iKataiainen et alj (l2007h suggest a magnetic moment // in the 
range 3 x 10^^ - 3 x 10^"^. The condition for accretion would then 
imply M > 5.3 x 10'^ g s"'. Although this value is probably at 
the high side, it might indicate that the a substantial contribution 
to the accretion luminosity comes from cyclotron radiation and 
not from the hard X-rays. 

Because of the wide range of tempe ratures found for the soft 
X-ray component. lAnzolin et al.l (l2008l) proposed that Tbb could 
be related to the magnetic field strength. Lower temperatures 
would indicate higher field systems because of larger irradiated 
WD areas by cyclotron radiation. However, the recent detection 
of sig nificant polarized emission in R XJ2133 dKataiainen et alj 
l2007ft and IRXS J173021.5-055933 (iButters et alJl2009h adds 
them to the subset of soft X-ray and polarized IPs together with 
V2400 Oph, PQ Gem and V405 Aur While the degree of polar- 
ization could also be affected by the accretion geometry, these 
five polarized systems possess soft blackbody components that 
span a wide range of temperatures, with three out of five display- 
ing hot blackbodies. The inferred magnetic field strengths are 
quite uncerta in for most systems, being 9-21 MG for PQ Gem 
(iPotter etalJ ll997). 9-2 7 MG for V2400 O ph (IVaethlll997h and 
~ 30 MG for V405 Aur jPiirola et al.ll2008l) . A 20 MG field was 



proposed for RXJ2133 ( IKataiainen et al.ir2007h . while no field 
estimate is given for IRXS J173021.5-055933 ("Butters et alj 
2009). To really test the proposal made by Anzolin et al. (20081, 
spectro-polarimetric measures of the magnetic field strengths are 
essential for all soft systems. Furthermore, there are indications 
that the reprocessed radiatio n emerges in the UV range al so in 
IPs, as found in PQ Gem ( Stavroviannop oulos et al.lll997h and 
UU Col (de Martino e t al. 2006.) . This aspect is essential to prop- 
erly determine the reprocessed energy budget. 



5. Conclusions 

The two CVs IGR0023 and RXJ2133 have been confirmed as 
true members of the IP class by using pointed X-ray observa- 
tions with XMM-Newton and Suzaku satellites. A strong pulsa- 
tion at the WD rotational period dominates the power spectra of 
IGR0023 below 2 keV, while for RXJ2133 it is detected up to 
12 keV. Both systems have fast rotating WDs, with P„ of 561 
s for IGR0023 and 571 s for RXJ2133. The latter is one of the 
most asynchronous systems among IPs, having a ratio between 
the spin and the orbital periods of 0.022. The fast rotation and 
very long orbital period, together with the detection of substan- 
tial polarized emission, suggests that RXJ2133 is likely a young 
IP that will evolve into a polar when attaining synchronism. 

Their broad-band X-ray spectra were analyzed including 
also INTEGRAL data, that allowed us to cover the energy range 
0.2-100 keV. These are well described by a multi-temperature 
plasma emission, with a minimum temperature of 0.2-0.3 keV 
and a maximum temperature of ~ 50 keV, which implies rel- 
atively massive WDs. While this might suggest that hard X-ray 
CVs harbor massive primaries, it has still to be understood if this 
is the only ingredient for a CV to be a hard X-ray source. 

In RXJ2133, a ~ 100 eV blackbody emission is also re- 
quired to fit the soft portion of the spectrum. The temper- 
ature of this component is similar to that found in V2400 
Qph ("de Ma rtino etal] \2QQ4) and IRXS J173021.5-055933 
(de Martino e t alj 1 2008 ). both of them also showing circularly 
polarized radiation. It is however different from that found in the 
other two soft X-ray and polarized IPs, PQ G em and V405 Aur 
dde Martino et alj|2004l: I Anzolin et ani2008l) . This casts doubts 
on the possible relation between the soft X-ray tempera t ure an d 
the magnetic field strength proposed by lAnzolin et al.l ( |2008|) . 
The present work also opens a further question on whether the 
reprocessed radiation in the soft and polarized IPs also emerges 
at UV wavelengths, as recently demonstrated by iKonig et akl 
(200&) for the polar prototype AM Her 

Furthermore, the fact that RXJ2133, V2400 Oph and 
IRXS J173021.5-055933 are bright hard X-ray sources and 
are also polarized IPs might not imply that cyclotron cool- 
ing decreases the hard X-ra y flux dWoelk & Beuermannlll996l : 
iFischer & Beu ermann 2001), at least for the magnetic strength 
values covered by these IPs. Observations aiming at measuring 
the magnetic field strength of IPs will help in making a clear 
picture of emission mechanisms in these systems. 
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Table 4. Spectral parameters of the best-fit models to the EPIC-pn, combined MOS and IBIS/ISGRI average spectra of IGR0023.The different models (see text) are 
indicated with A, B and C in the first column. Errors indicate the 90% confidence level of the corresponding parameter. 





NhI 










(keV) 




kT2 

(keV) 




kT3^ 
(kev~ 




E.W. ■' 
(eV)" 


xi 

(^2/d.o.f.) 


A 
B 
C 


1.79: 

2.0 ± 
2.0 ± 


±0.07 
0.1 
0.1 


0.8 ±0.2 


0.44 ± 0.04 
0.36 ± 0.06 


0.7 ±0.2 
0.5 ±0.1 

' -0.2 


17+"°^ 

^ ' -0.04 

17+"°^ 


3-0!?:? 


14 ±2 


•'■^-0.3 
'^■'--0.9 


28!|^ 
> 27 


^■'-0.5 

0.3 ±0.1 


100+^" 
100!g 

94+24 


0.86 (892/1033) 
0.86 (886/1032) 
0.85 (870/1030) 



" Column density of the total absorber in units of 10^' cm"-. 
' Column density of the partial absorber in units of 10^^^ cm"^. 

Covering fraction of the partial absorber. 

Metal abundance. 

Index of the power-law emissivity function of the CEMEKL component. 
^ Normalization constant of the first optically thin plasma component in units of lO""*. 

* Normalization constant of the second optically thin plasma component in units of 10"\ 

* In model A this the maximum temperature of the CEMEKL. 

' Normalization constant of the third optically thin plasma component (CEMEKL in model A) in units of 10"^. 
' Equivalent width of the 6.4 keV Fe line. 



Table 5. Spectral parameters of the best-fit models to the EPIC-pn, combined MOS and IBIS/ISGRI average spectra of RXJ2133 for May 2005 (7) and July 2005 
(2), as well as to the XIS and HXD spectra (i). The different models (see text) are indicated with A, B and C in the first column. Errors indicate the 90% confidence 
level of the corresponding parameter. 



E.W. * 
(eV)" 



A'hI 



7+m 

-0.2 

2 1.7 (fixed) 

3 1.7 (fixed) 



56+"-"^' 

"•-'"-0.04 

0.58+!!!;3 



(eV) 



"94T 
97+f 
102 ± 1 



6+"-^ 

"•"-0.1 
"•"-0.1 



(keV) 



kT2 
(keV) 



(keV) 



55^ 

45+23 



(X- 1 d.o.f ) 



1 1.7+ 



1 3+"-^ 
^•■-'-0.1 

1.1 ±0.2 



1.06 ± 0.07 0.48 



7+0.12 

-0.08 

2 1.7 (fixed) 

3 1.7 (fixed) 



-0.02 



4.9^ 
3.7 



0.5 



:0.2 0.8 ±0.1 



■)+o:u4 

-0.02 

0.59 ± 0.04 

£;+0.02 
-0 01 



^•^-0.6 

4 3+"-5 

^•-'-0.4 

3.1 ±0.1 



9+" - 

"•^-0.1 

0.9 ±0.2 

1 (fixed) 



3.5: 

3.5 : 



:0.4 

0.6 



1+0.4 

-0.2 

4.6!" 
7 3+^ '^ 

-^•^1.4 



5+03- 

-0.3 

5.1+'-' 
3 6+'^-^ 

^- '-1.7 
1.0 



"27T 

> 71 

"395^ 

39!^^ 

> 61 



■^•°"-o.08 



159 
149^ 
181 



±22 

16 
31 

± 10 



1 70+"^"^ 

^•'"-0.07 

1.3 ±0.2 

1.46+!;?! 



0.98 (1416/1444) 
1.05 (819/780) 
1.12(2907/2608) 



140 ± 30 
173 ±9 
1565^ 
142+^' 
170 ±9 



1 1.7+ 



1.4 ±0.2 

1.2 ±0.2 

1.03 ± 0.07 



0.52+ 



0.46+ 



98 ± 
103 : 
107 : 



0.5 ±0.1 

0.7+" 

0.9+?? 



1.4+ 



0.9+ 



96+^ 
99 + 2 
100+^ 



5 4+U.« 

''■^-0.9 
4.7 ± 0.5 

fl.O 

-0.8 



1 +0.3 

"-0.2 
q+0.2 
^-0.1 



0.99(1431/1443) 
1.06 (822/779) 
1.08 (2818/2606) 



1 1.7 ±0.1 

2 1.7 (fixed) 

3 1.7 (fixed) 



1.3 ± 
1.1 ± 
1.10 



0.2 
0.2 
t0.09 



0.54: 
0.57: 



0.03 
0.04 



45+0'" 

"• -^-0.02 



4.7+ 



0.7 ± 0.2 

"• '-0.2 
9+°-2 

"•^-0.1 



4 1+ + ' 

-1-7 

1 6+"-' 

-06 



6.6 ±1.0 3.5 



1.46!;» 

l^3!!Ji 

1 50+""^ 



0.97 (1403/1441) 
1.04(805/777) 
1.07 (2794/2603) 



" Column density of the total absorber in units of 10^' cm"^. 
' Column density of the partial absorber in units of 10^^^ cm"^. 
Covering fraction of the partial absorber. 

Normalization constant of the BBODY component in units of lO"'*. 
Metal abundance. 

^ Index of power-law emissivity function of the CEMEKL component. 

* Normalization constant of the first optically thin plasma component in units of lO"''. 

'' Normalization constant of the second optically thin plasma component in units of 10"^^. 
' In model A this the maximum temperature of the CEMEKL. 

' Normalization constant of the third optically thin plasma component (CEMEKL in model A) in units of 10"^. 

* Equivalent width of the 6.4 keV Fe line. 



